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abstract 

Rose, Ferrante. and Smith (ref. 7) have discovered scaling relations which map 
the adhesive binding energy of Ferrante and Smith (ref. 6) onto a single universal 
binding energy curve. The energies in ref. 6 are calculated for all combinations of 
Al(lll), Zn(OOOl). Mg(OOOl), and Na(llU) in contact. The scaling involves 
normalizing the energy to the maximum binding energy and normalizing distances by a 
suitable combination of Thomas-Fermi screening lengths. A simple mathematical 
expression is found to accurately represent the universal curve. E* (a*) - -(1 + ea*) 
exp (>Ba*) where E* is the normalized binding energy, a* is the normalized 
separation, and a is the fitting parameter. Rose et a1. (ref. 7) have also found 
that the calculated cohesive energies of K. Ba, Cu, Mo, and Sm scale by similar 
simple relations suggesting the universal relation may be more general than for the 
simple free electron metals for which it was derived. In this paper we outline these 
results and discuss them in relation to topics cf interest in adhesion, friction, and 
wear. 


INTRODUCTION 

An important aspect of friction and wear of metals is the nature of the adhesive 
force between the surfaces. Bowden and Tabor (ref. 1) have used adhesion between 
metals to explain parts of the friction and wear process. Of particular interest in 
these processes is the physical nature of the forces in dissimilar metal contacts. 
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There have been numerous attempts to postulate the nature of these forces but few 
model calculations. For example, adhesion between d1:slm11ar metals was explained In 
terms of solubility (ref. 2) but without accompanying mocel calculations. Buckley 
(ref. 3) has shown that metals with low mutual solubility can have higher adhesion 
than those with high solubility. 

This lack of model calculations has occurred because of the theoretical diffi- 
culties Involved In treating solid surfaces. In the past decade the electron density 
functional formalism of Hohenberg, Kohn and Sham (ref. 4) have been used to success- 
fully calculate the surface energies (ref. S) of simple metals. Ferrante and Smith 
(ref. 6) recently applied this formalism to the problem of metallic adhesion for both 
similar and dissimilar metals In contact. Rose, Ferrante and Smith (ref. 7) found 
that scaling laws exist which map these results onto a universal adhesive energy 
curve. In this paper, we discuss the results of these calculations and their Impli- 
cations to a more general category of materials, the transition metals, and then sum- 
marize the current state of adhesion theory. 


USl OF KOHN-SHAM FORMALISM FOR ADHESION CALCULATIONS 


The calculational formalism and methods used for obtaining self-consistent interface 
electronic structure will now be presented. A much more extensive description is 
given In refs. 6 and 8, particularly of numerical techniques. 

The adhesive interaction energy, E^^, between two metal surfaces is a function 
of the distance between the two surfaces, a (see Fig. 1). ^ad defined as the 
negative of the amount of work necessary to increase the separation from a to « 
divided by twice the cross-sectional area A. Thus 

- £(«>]/2A, U) 

where E is the total energy. For identical metals, E^^ is the negative of the 
surface energy when a is at the energy minimum. 

According to the density functional formalism of Hohenberg, Kohn, and Sham, 

(ref. 4) the total energy is given by (atomic units are used throughout unless other- 
wise specified) 
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v(r) is the ionic potentle), and ni^) is the electron number density. The first two 

terms in £q. (2) are the electron-ion and ion-ion interaction energies* respec- 

tively. z is the ionic charge and is the distance between ion core nuclei 
(there is no ion core overlap in the systems considered here). T^(n(D} is the 
kinetic energy of a system of noninteracting electrons with the same density n(»^, 
the next term is the classical electron-electron interaction energy, and is 
the exchange-correlation energy. 

for metals like 2n, Ng. Al, and Na» the Jellium model (Fig. 1) is a good 
zeroth-order approximation, and the difference between the total pseudopotential and 
the potential due to the jellium is small for the closest packed plane. Thus for a 

given separation a, one obtains E to a first-order perturbation approximation as 


Ein(r)} « A J Vj(y,a)n(y,a)dy ♦ j ^ ♦ F{n(y.a)> 

i^j 

+ A / av (y.a)n(y.a)dy (4) 

where Vj is the potential produced by the jellium, y is the direction normal to 
the surfaces, and av is the average, over planes pa*allel to the surface, of the 
difference in potential between an array of pseudopotentials and the jellium. 
Following Lang and Kohn (ref. b) the Ashcroft pseudopotential is used: 
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where is determined empirically and is close to the ion core radius. The 
electron number density is obtained from a set of self-consistent equations. 


The bchroedinger equation: 
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and Mith Poisson's equation 

. -4i[n(y,a) - n^(y)] (^) 

d/ 

^(y«a) Is t«)e electrostatic pot -it a) n^(y) Is the jelllun density (Fig. 1} 
kp Is thc: Fertnl wave vector magnitude and t|t^^^(y) Is the wave function 
(1 - L, R; I.e., left and right). The calculation Is broken up Into two regions, 
below the conduction band of the less dense metal, and above 11 where the solutions 
are doubly degenerate. 

It Is useful to combine the first two terms of Eq. (4) along with the classical 
electron-electron Interaction term of F{n(y,a)} as follows: 


- j A y* p{y,a)p(y,a)dy ♦ (8) 

where e(y,a) Is the net charge density of the zeroth-order jelllum solution. 

Is the exact difference between the Ion-Ion and the jelllum-jelllum inter- 
action. It Is shown In ref. 9 that Is negligible unless the facing planes 

are in registry, I.e., commensurate. In this calculation, we assume Incommensurate 
adhesion - 0), since registry Is not obtained with dissimilar metals In con- 

tact without corresponding loss of energy due to strains In the lattice which would 
be difficult to evaluate. 

The exchange-correlation energy E^. Is written In the local-density approx- 
Imatlon (LOA), 


• J' n(r)e^^(n{r))dr ^3) 

where Cj^^(n(r)) Is the exchange-correlation energy of a uniform electron gas of 
number density n(r). We use Wigner's Interpolation formula for the correlation 
energy and the Kohn-Sham exchange energy. 

It remains to specify the kinetic-energy functional of Eq. (3): 
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where 


t5{n(y,a)} 


t • 'Si 


(occ.) 
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♦ [vgf^(n;H - Vg^|;(n;y)]n(y,a) (11) 

The sum 1$ over all occupied states. The sunmatlon index 1 again refers to 
degenerate wave functions as In Eq. (6). In ref. 6. an expression for 
T (n(y,a)} - T^{n(y,U)} Is presented which Is based on kinetic-energy 
densities In the Interface region. We find this approach more natural for our prob- 
lem, as It Is the Interface region In which the large changes In kinetic energy den- 
sity occur upon adhes1o^. 

The results of this calculation for all dissimilar combinations of Al(lll), 
ii..(0001), Mg(OOOl) and Na(llO) are shown In Fig. Z. as can be seen there Is a wide 
variation of shapes and binding energies and the curves look very similar to binding 
energy curves for diatomic molecules. 

Scaling rules for adhesive energies 

The calculations needed to obtain the results of Fig. i are quite difficult and 
time-consuming. It Is of interest, therefore, to look for some general similarities 
between these curves which may generalize the results. Hose, et al. (ref. 7) have 
found such scaling laws. It will now be shown that the curves of Fig. z, as well as 
those of the Identical metal contacts (ref. 6) (A1(1U) - A1(U1), etc.), can be 
simply scaled Into a universal curve. This scaling is motivated by the expectation 
that metals having shorter screening lengths would have adhesive energy curves which 
rise faster with separation. That Is, the metals would screen the disturbances 
caused by creating the surface over a shorter distance. This suggests that for Iden- 
tical metal contacts the separation Is scaled by the Thomas Fermi screening length 
X - (9n/4)^^'^ r^ au, where the bulk electron density n^ ■ 3/4irr^. 

When we encounter bimetallic contacts as represented in Fig. 1, a length scaling 
appropriate to both metals must be considered. In that case, we chose to scale by an 
arithmetic average, The energy amplitude was scaled by Its 

equilibrium value, aE = where a^ Is the equilibrium separation. 

Explicitly 

Here Is the univeral adhesive energy function and a* • 2(a - 

Figure 3 shows the results of scaling the calculated adhesive energies. An 
analytical fit Is given by • - (1 ♦ aa*) exp (-aa*) with a • 0.90. The 



universollty of tho scaled adhesive energy corve Is truly remarkable. One can see 
that the results for all ten bimetallic contacts lie very close to the universal 
curve. This Is true even though the bulk metallic densities In the various metals 
vary by a factor of eight. 

It Is Important to understand how the fortunate result of a universal energy 
curve comes about. In the following, we attempt to provide a plausibility argument 
within the jelHum model. First, we have found that sol1d»vacuum density d1s> 
tributlons, n(y) (ref. 10), scale rather accurately with x. That Is, there Is, to a 
good approximation, a universal number density distribution n*(y - a/2) where 

n(y - a/2) - n^n*[x“^y - a/2)] (13) 

Here a/2 1$ the coordinate of the jelllum surface. There a similar scaling for 
the Kohn>Sham effective one electron potential: 

Vf^^ " • ''b 

We were motivated to look for this scaling by the fact that the Thomas Fermi equation 
scales (refs. 11,12) exactly with y In units of x. Secondly, we have found (ref. 
6) that the total numbi^r density in the bimetallic interface is given to a fair - 
accuracy by a simple overlap of the corresponding solid-vacuum distributions. This, 
and the stationary property of E[n] indicates that it would be a good approximation 
to use overlapping solid-vacuum number density and potential distributions. 

Thus, in a first order perturbation approximation, we have for identical metal 
contacts: 

^00 

Combining Eqs. (3) to (5), we have: 

tad^*J ■ (l/«)(n^Vg) y* n*(y - a*/2)v*^^(y ♦ a*/2)dy (16) 

The integrand in Eq. (16) is independent of r^. The constants in front of the 
integral are independent of a, and thus Eq. (16) scales exactly as we scaled the 
adhesive energy curve to give Fig. 3. Although our plausibility argument is 
restricted to Jelllum interfaces, we note that the calculated adhesive energies 
include the ion-ion term exactly for a rigid lattice model and the electron-ion term 
in first order perturbation theory. 

UiSCUSSlON 

There are several topics of importance in metallic adhesion: the nature of the 
attractive forces, the range of such forces, and finally, are there any 
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generillzdtions that can be found t The work of Ferrante and Smith (ref. 6) indicates 
that electron sharing similar to covalent bonding In molecules 1:. sufficient to give 
strong bonding. They also determine the magnitudes of such forces for several simple 
metals (table 1). These results also establish that an approximate range for the 
forces Is an Interplanar spacing. Finally, Rose, Ferrante and Smith (ref. 7) have 
found scaling laws which generalize these results. 

In metallic adhesion of real surfaces further considerations come to mind. 

First, real surfaces are not flat and consequently, how Important are the strong 
forces compared to the weaker, long-range van der Uaals forces? Next, the results of 
these theoretical calculations represent ideai strengths. Real contacts fall at 
conditions lower than at Ideal strengths. Finally, metals used In most engineering 
applications are either transition metals or alloys. 

inglesfleld (ref. 13) examined the problem of van der Waals vs. strong Inter- 
action forces using the calculations of Ferrante and Smith (ref. 9). He concluded 
that In spite of the fact that the strong forces only dominate at the positions of 
asperity contact, these strong Interaction forces dominate the strength of the 
contact. 

The fact that the strength of real contacts Is dominated by defect structures 
and the mechanical properties of the solids is an issue that ultimately must be 
addressee. Also, these results apply for brittle fracture; however, metals will 
undergo ductile extentlon before fracture. Since adhesion and other contact pro- 
cesses are complicated combinations of physical and mechanical properties, each part 
of the process must be understood Independently. This study has been car. led out In 
this spirit. The results of table 1, however, contain some aspects of these con- 
siderations. A question of interest is whether a dissimilar contact is weaker, 
stronger or in-between contacts between the same metal. The first column In table 1 
(labeled Wint ■ 0) gives the binding energy of a contact when perfect registry Is not 
obtained. Perfect registry can only be obtained when the same crystallographic 
planes of the same metal are in contact. A more complete discussion of this issue Is 
given in ref. 9. The second column gives the binding energy in the ideal case of 
perfect matching of planes which doesn't occur in a real contact. Examining table 1, 
it can be seen that the Interfacial energy Is slightly lower than or comparable to 
the bulk energy except for the case of sodium, where the interface between sodium and 
any of the other metals is clearly stronger than the bulk. In a real contact, the 
two surfaces would probably distort to attempt to come into registry with the for- 
mation of defect chains to accomodate the distortions, thus modifying the Interfacial 
energy. In ge'^eral, the interfacial binding energy is comparable to the energy of 
the weaker partner In tht different metal contacts examined. 

A final concern with these results is the limitations of the calculations. The 
quasi-three-dimensionality of the adhesion calculations limits the reliability of the 
calculations to the densest packed planes of the simple metals. 
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Most tnglneerli.g naterlrU are transition Mtals. To properly handle transition 
metals theoretically would require extensive three-dimensional calculations which 
would be quite difficult at this stage. Yanlv (ref. 14) and Allan* Lannoo and 
Dobrzhinski (ref. IS) have applied the tight-binding approximation to the transition 
metal problem. Although these calculations probably give accurate trends, the 
approximations used make quantitative results questionable. 

The question arises concerning whether the results are more general than for the 
simple metals examined. Recently, theoretical binding energy curves have become 
available for several bulk metals [Carlsson et al. (ref. 16) (Mo, K, and Cu) and 
Herbst (ref. 17) (Sm Sro and Ba)]. These total cohesive energy curves were 
calculated as a function of the separation between atoms for a uniformly dilated 
lattice. We characterize the density of the lattice In terms of the Wigner-Seltz 
radius, r^^ ■ (3/4im^)^^^ where n^ Is the atom density. As shown In fig. 

4, these quite disparate cohesive energy curves can be approximately scaled Into a 
universal function, £^, which Is also defined In Eq. (12) If we replace 
by everywhere. a£ Is the cohesive energy at the equilibrium spacing r^^^^ 
and a* • (r^^ - where x Is again the Thomas-Fermi screening length. 

The value of r^ used to determine x was determined by the equilibrium inter- 
stitial electron density (refs. 17,18). The binding energy of Ho, K, Ba, 

Sm [4f (5d,6s)^] and Sm*^[4f^(5d,6s)^] fall closely on a single curve 
with a a 1.16 where we have used the same analytic form as for the adhesive ener- 
'iles. The value of a differs from that appropriate for adhesive energies presum- 
ably In part because all atoms change their positions in the bulk cohesive energy 
calculations while the adhesive energy curves assume that atomic planes are moved 
rigidly. The results for Lu has the same shape, but a somewhat different a than 
the other metals. We do not understand this variation. 

The cohesive energy calculations of Carlsson et al. (ref. lb) (ASM density func- 
tional theory) and Herbst (ref. 18) (relativistic Hartree-Fock ) are quite different 
from each other and from the perturbative-density functional results of ref. 6 for 
^ad* nature of cohesive bonding In these metals Is quite varied. Ba is a 

divalent band overlap metal; Sm Is an f-electron metal; Ho and Cu have Important 
d-band Interactions; while K Is a simple metal. That such different metals cal- 
culated In quite different ways fall on a single curve Indicates the generality of 
the scaling relations. We note that the "tail" region of the screening charge dis- 
tribution around metal Ion cores can be represented by electron gas parameters. It 
Is known that the screening charge density and potential In this region scale as 
r^^^ to a fair approximation (ref. 19). Then the argument of £qs. (13 to 16} 
applies to the long-range Interaction of screened bulk Ions with i a/2 referring 
to the relative lattice positions and y replaced by r*. 
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The use of the electron density at is reminiscent of the approach of 
Niedema et ai. (ref. 20) whose work complements our own. They calculate equilibrium 
energies (e.g., heats of formation). We. on the other hand, use similar equilibrium 
quantities to determine the form of the binding energy-distance relation. 

The plotted energies of Fig. 4 were computed for the same electronic con- 
figuration at all atomic separations as were the interface calculations. The 
adhesive energy scaling was illustrated for simple metals. However, the appearance 
of an analogous scaling relation for the bulk energy of transition metals indicates 
that the adhesive energy scaling may well extend beyond simple metals. 

In conclusion we have discovered scaling relations which map adhesive energies 
onto a universal binding energy curve. Similar scalings for cohesive energies for a 
wide range of metals suggest that this relation may apply to more than Just the 
simple metals. Finally, a simple analytic energy-separation scaling relation is pro- 
vided which may be of use in practical analysis of contacts. 
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Netal combination Binding energy 
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ergs/cro* registry 


Ai(iii)-Ai(in) 

490 

715 
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505 
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Mg(001)-Hg(0001) 

Na(U0)>Na(110) 

460 
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A1(M1)-Zn(0001) 

520 


Ai:Ul)-Mg(0001 

505 


A1(ill)-Ka(110) 
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Zn(0001)-Mg(0001) 

490 


Zn(0001)-Na(110) 

325 


Mg(0001)-Na(110) 

310 





ADHESIVE ENERGY, ergfcm 



POSITION y. nm 


Figure 1. - Electron numtier 
density n and jelllum Ion 
charge density n^ for an 
Al-Mg contact. 







Figure 2. • Adhesive binding energy versus the separation a between 
th ) surfaces Indicated. 




SCALED SEPARATION, a« 

Figure 3. - Adhesive energy results from Fig. 1 above and Fig. 4 of 
Ref. 1 scaled as described tn the tex^ a« - 2(a - * A2^* 



SCALED CHANGE IN WIGNER SEITZ RADIUS. a« 


Figure 4. - Bulk energy of various metals scaled and described 
In the text: a* - (r^ - r^^MX. A value of p • 1. 16 was 
obtained by a least squares fit to the data. 




